Extracting the condensate density from projection experiments with Fermi gases 
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A debated issue in the physics of the BCS-BEC crossover with trapped Fermi atoms is to identify 
characteristic properties of the superfluid phase. Recently, a condensate fraction was measured 
on the BCS side of the crossover by sweeping the system in a fast (nonadiabatic) way from the 
BCS to the BEC sides, thus "projecting" the initial many-body state onto a molecular condensate. 
We analyze here the theoretical implications of these projection experiments, by identifying the 
appropriate quantum-mechanical operator associated with the measured quantities and relating 
them to the many-body correlations occurring in the BCS-BEC crossover. Calculations are presented 
over wide temperature and coupling ranges, by including pairing fluctuations on top of mean field. 

PACS numbers: 03.75.Hh,03.75.Ss 
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The current experimental advances with trapped Fermi 
atoms have attracted much interest in the physics of the 
BCS-BEC crossover. In this context, one of the most 
debated issues is the unambiguous detection of super- 
fluid properties on the BCS side of the crossover. Several 
attempts have been made in this direction. They in- 
clude absorption images of the "projected" density pro- 
files for ""^K and ^Li 0, rf spectroscopy to detect 
single-particle excitations j^l, and measurements of col- 
lective modes 0, 13- In addition, a number of schemes 
to detect superfluid properties on the BCS side of the 
crossover_ have been proposed, including Josephson oscil- 
lations 6] and vortices Y]. 

In particular, the experimental procedure of Refs. 0,0 
pairwise "projects" fermionic atoms onto molecules, by 
preparing the system of trapped Fermi atoms on the BCS 
side vifith a tunable Fano-Feshbach resonance and then 
rapidly sweeping the magnetic field to the BEC side. In 
this way, the same two-component fit of density profiles 
routinely used for Bose gases is exploited to extract from 
these "projected" density profiles the analog of a conden- 
sate fraction, which is now associated with the equilib- 
rium state on the BCS side before the sweep took place. 

Purpose of this paper is to provide a theoretical inter- 
pretation of the experiments of Refs. P, 0] i by obtaining 
the "projected" density profiles in terms of the correla- 
tion functions of the Fermi gas at equilibrium. This will 
be based on a number of physical assumptions which we 
associate with the experimental procedure of Refs. 
Our calculation evidences how the projection procedure 
amplifies the visibility of the emergence of the condensate 
as the temperature is lowered below the critical temper- 
ature Tc, when compared with the ordinary density pro- 
files of Ref . . We also attempt an analysis of the "pro- 
jected" density profiles in terms of a two-component fit, 
in analogy to what is done with the experimental data 
Q, Q . A prediction is further made of a reduced molecu- 
lar fraction that depends on the initial many-body state, 
in agreement with a late experimental evidence [2| . 

Inclusion of pairing fluctuations on top of mean field 



along the lines of Ref. JS^ enables us to cover a wide 
temperature range even in the intermediate- and strong- 
coupling regimes, in contrast to Refs. 0,0] where only 
mean field was taken into account. 

To account for the "projected" density profiles of 
Refs. 0,I3> we consider the boson- like field operator: 

vI^B(r)= jdpcjy{p)i:^{v-p/2)i,^{v + p/2). (1) 

Here, is a fermion field operator with spin a, and 

(j){p) a real and normalized function which specifies the 
probability amplitude for the fermion pair. An operator 
of the form |^ was considered in Ref. 0| to obtain the 
condensate density for composite bosons. 

Our theoretical analysis of the experiments of Refs. P, 
01 is based on the following physical assumptions, that 
we infer from the experimental procedure: 

(i) Atoms of a specific spin state were detected, which 
originate from the dissociation of molecules after apply- 
ing an rf pulse. The object of the measurement is thus 
the bosonic (molecular) density ^^(r) at position r (and 
not the fermionic (atomic) density n(r)). 

(ii) A rather large conversion efficiency into molecules re- 
sults when rapidly sweeping the magnetic field in the ex- 
periments. This suggests that molecules form just past 
the unitarity limit on the BEC side, where the "final" 
molecular wave function and the many-body correlations 
for the "initial" states considered in Refs. ffl,^ have maxi- 
mum overlap. Correspondingly, we assume that the wave 
function in the expression ^ refers to this "final" cou- 
pling, and represent it by 4'fip)- As molecules form in a 
medium, we take into account the effect of Pauli blocking 
in analogy with the original Cooper argument 0| and 
identify (f> f (p) with the bound-state solution of the two- 
body problem, with the condition that its Fourier trans- 
form (j)f(M) vanishes when the magnitude of the wave vec- 
tor k is smaller than a characteristic value k^^ = ^2mpf 
for nf > ( m being the fermion mass), while no con- 
straint is enforced for pf < 0. Here, the value of the 
chemical potential pf depends on the "final" coupling 
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at which the molecular state is assumed to form, thus 
implying that some sort of local equilibrium can be es- 
tablished around a molecule. We shall present our calcu- 
lations for two "final" couplings that bound the interval 
where the maximum overlap occurs, and are at the same 
time representative of the two cases where /i/ is positive 
or negative. 

(iii) In the experiments, bosonic Thomas- Fermi (TF) pro- 
files for the molecular condensate were extracted from 
position-dependent density profiles, thus entailing an as- 
sumption of thermal equilibrium. We assume that this 
thermal equilibrium corresponds to the state prepared 
before the rapid sweep of the magnetic field. The valid- 
ity of this assumption is supported by a recent exper- 
imental stiidy of the formation time of a fermion-pair 
condensate|9|. In our calculations, we then use (• • ■)i as 
expressions for the thermal averages, where the suffix i 
stands for "initial" . 

All these assumptions are summarized by stating that 
the "projected" bosonic density profile given by 

ni^ir) = (vl//(r)tvl,/(r)}, (2) 

represents the "in situ" molecular density which would 
be measured after the rapid sweep but before the cloud 
expansion performed in Refs. P, Q (connection with the 
density measured after the expansion will be discussed 
below). In this expression, the boson-like field opera- 
tor of Eq. ^ contains the final molecular wave function 
4>f{p) on the BEC side of the crossover, while the ther- 
mal average {■ ■ •)i is taken with reference to the state in 
which the system was initially prepared. 

Consistently with our previous work Til, we describe 
the interaction term of the many-fermion Hamiltonian 
via an effective single-channel model. The parame- 
ter {kpaF)~^ then drives the crossover from the BCS 
side (identified by {kpap)~^ ^ —1) to the BEC side 
(identified by 1 ^ (kpaF)"^) across the unitarity limit 
(kpap)~^ = 0. Here, ap is the two-fermion scattering 
length and the Fermi wave vector kp results by setting 
kp/{2m) equal to the nonintcracting Fermi energy. 

The calculation proceeds by expressing the four- 
fermion field operator in Eq. |5J| in terms of 
the two-particle Green's function ^72(1, 2, 1', 2') = 
(r^[*(l) ^-(2) *t(2') *t(i')])^ where is the imaginary- 
time ordering operator. We have introduced the spinor 
\I'(r) = (■0|(r), '0|(r)) as well as the short-hand notation 
1 = (ri , Ti , £i ) with imaginary time r and spinor compo- 
nent £, such that ^'(1) = exp{7^Ti}*£j (n) exp{-i^Ti}. 
The thermal average contains the grand-canonical Hamil- 
tonian K = H — fiN with fermionic chemical potential fi 
and is taken in the initial state, as specified above. 

The two-particle Green's function Q2 is, in turn, ex- 
pressed in terms of the many-particle T-matrix, by solv- 
ing formally the Bethe-Salpeter equation as follows: 

^2(1, 2, 1', 2') = g{i, 1') g{2, 2') - g{i, 2') g{2, 1') 



- y"d3456 ^(1, 3) g(6, 1') r(3, 5; 6, 4) ^(4, 2') 0(2, 5) (3) 

where ^(1,1') = - (T^[*(l) ^-+(1')]) is the fermionic 
single-particle Green's function. Accordingly, the "pro- 
jected" bosonic density Q reads: 

ni\r) = J dpdp' 4>f{p) cpfip') g^(l,2,l',2') (4) 

where 1 = (r- p/2, r + 277, ^ = 2), 2 = {r + p' /2,t,1= 1), 
V ^[v + p/2,T + ir],£^l), and 2' = (r - p72, T + 77, ^ = 
2) (?7 being a positive infinitesimal). 

Implementation of the above expressions to the 
trapped case is readily obtained via a local-density ap- 
proximation, whereby a local gap parameter A(r) is in- 
troduced and the chemical potential /i is replaced (when- 
ever it occurs for both "initial" and "final" couplings) by 
the quantity /i(r) = /i — V{t) that accounts for the trap- 
ping potential V{y). 

The three terms on the right-hand side of Eq. Q corre- 
spond to physically different contributions to the expres- 
sion Q. In particular, the first term can be written as 

ja/ip, where a/i — Jdp(f)f{p)g[2 {p,t = —v) represents 
the overlap between the fermionic correlations (embodied 
in the anomalous single-particle Green's function ) 
and the molecular wave function (pf. This contribution 
vanishes with the gap parameter A when approaching Tc, 
and is identified with the condensate density for compos- 
ite bosons when both the "initial" thermal equilibrium 
and the "final" molecular wave function are taken at the 
same coupling deep in the BEC region Only this 

contribution was considered in Ref. in connection 
with the experiments of Refs. 

The second term on the right-hand side of Eq. ^ rep- 
resents fermionic correlations in the normal state, that 
are relevant in the presence of an underlying Fermi sur- 
face. This term is most sensitive to the "final" molecular 
wave function (pf being affected by Pauli blocking when 
/i/ is positive. This term would be irrelevant if the "ini- 
tial" thermal equilibrium and the "final" molecular wave 
function were taken deep in the BEC side. When /iy 
is positive, this term can lead to an overestimate of the 
value of the condensate fraction when the "projected" 
density profiles are fitted in terms of TF and Gaussian 
functions, as argued below. Both this and the previ- 
ous contribution were considered in Ref. (where the 
"final" coupling was, however, taken deep in the BEC 
region) . 

The third term on the right-hand side of Eq. J^Jl will be 
calculated in the following within the off-diagonal BCS- 
RPA approximation considered in Ref. • This contri- 
bution is identified with the noncondensate density for 
composite bosons when both the "initial" thermal equi- 
librium and the "final" molecular wave function are taken 
at the same coupling deep in the BEC side [isj . It is thus 
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FIG. 1: 



Axially-integrated "projected" density profiles (in 
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units of Rp wfiere _Rf is tlie fermionic Tliomas-Fermi ra- 
dius) for three "initial" coupling values and for two "final" 
couplings: 0.40 (upper panel) and 1.50 (lower panel). 
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FIG. 2; Ratio Nmoi/ {Nf /"i) vs {kpap)^^ for two tempera- 
tures and "final" couplings. 



of particular importance for increasing temperature when 
approaching the normal phase. 

When i = f deep in the BEC region, the "projected" 
density profile ?^^(r) coincides with (half) the ordinary 
density profile n(r) calculated at the same coupling, fn 
the following, the values for the local chemical potential 
and gap parameter, to be inserted in the expression Q 
for n;^'(r), are obtained with the theory of Ref. '§] where 
pairing fluctuations are included on top of mean field. 

Figure 1 shows the axially-integrated "projected" 
density profiles calculated for the coupling values 
{kpaF)'^ = (-0.50,-0.25,0.00) and for the two repre- 
sentative values 0.40 (upper panel) and 1.50 (lower panel) 
of the "final" coupling [kpap)^^ ■ Two characteristic 
temperatures (just above the critical temperature and 
near zero temperature) are considered in each case. Note 
the marked temperature dependence of the "projected" 
density profiles when entering the superfiuid phase, as 
signaled by the emergence of a "condensate" component 
near the center of the trap. This contrasts the milder de- 
pendence (especially on the BCS side) of the density pro- 
files without projection 0. The "projection" technique 
introduced in Ref. jlj is thus demonstrated to amplify 
the effects due to the presence of a condensate on the 
density profiles, which would otherwise be almost tem- 
perature independent on the BCS side of the crossover. 

In Fig. 1 the densities are normalized to half the to- 
tal number Np of fermionic atoms. This number dif- 
fers, in general, from the total number Nmoi of molecules 
obtained by integrating the "projected" density profiles. 
In particular, Nmoi can vary significantly when scanning 
the "initial" coupfing {kFaF)i^ on the BCS side of the 
crossover for given "final" molecular-like state. This ef- 
fect is shown in Fig. 2 for the same temperatures and 
"final" couplings of Fig. 1. Our finding that the total 
number Nmoi of molecules constitutes only a fraction of 
the original atom number Nf/2 for each spin state is 
supported by the experimental results of Refs. 0] and 0- 
In addition, our prediction that the reduced value of the 



molecular fraction depends on the "initial" many-body 
state is in agreement with the experimental evidence re- 
cently reported in Ref. |^. Note that, for both values 
of the "final coupling", the total number of molecules 
increases upon lowering the temperature. This result in- 
dicates that the conversion efficiency for the condensate 
fraction is larger than for the thermal component, as also 
observed in Ref. . 

In our procedure, the condensate and noncondensate 
components of the "projected" density profiles are cal- 
culated separately. By our definition, they correspond 
to the first term and to the remaining terms on the 
right-hand side of Eq. (j^J, respectively. The conden- 
sate fraction is obtained accordingly from the ratio of 
the corresponding areas. Yet, the total "projected" den- 
sity profiles obtained theoretically could be analyzed in 
terms of a two-component fit with a TF plus a Gaussian 
function (or, better, a g3/2 function for the Bose gas), 
in analogy to a standard experimental procedure. This 
kind of analysis is reported in Fig. 3 for the two low- 
temperature curves shown on the right panels of Fig. 1. 
In both cases, a good overall fit is obtained by the 
square method. Separate comparison is also made in the 
figure with the theoretical condensate and noncondensate 
components defined above, which appears rather good for 
the value 1.50 of the "final" coupling while for the value 
0.40 an overestimate (of about 50%) of the condensate 
results from the fit. This discrepancy stems mostly from 
the second term on the right-hand side of Eq. Q , which 
contributes to the TF component of the fit owing to a 
peculiar shape of the corresponding "projected" density 
profile. 

Extention of this analysis to "initial" couplings on the 
BCS side reveals unconventional forms of the theoretical 
condensed and noncondensed contributions to the "pro- 
jected" density profiles, so that the above two-component 
fit fails. For negative values of the initial coupling when 
the two-component fit fails, we have verified that the dif- 
ference Nmoi{T = 0) — Nmoi{T = l.lOTc) approximately 
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FIG. 3: Two-component (TF plus g3/2) fit of the axially- 
integrated "projected" density profiles (in units of i?^^) at 
the unitarity hmit for T/T^ = 0.05. The TF (dashed line) 
and g3/2 (dotted line) components of the fit are compared 
with the theoretical condensate (full line) and noncondensate 
(dashed-dotted line) components. 



coincides with Nq{T = 0) within a relative error not 
larger than 10% when {kpap)^^ = 0.4. This suggests 
a practical prescription to extract Nq{T = 0) from the 
values of Nmoi at low temperature and slightly above the 
critical temperature, without relying on a two-component 
fit. 

In this respect, recall that our theoretical "projected" 
density profiles are calculated when molecules just form 
on the BEC side near the unitarity limit. In the exper- 
iments, however, a further ramp of the magnetic field 
is performed together with a subsequent cloud expan- 
sion. Only at this stage the profile of the cloud is de- 
tected. Comparison between theory and experiments 
is thus meaningful since the further ramp of the mag- 
netic field and the subsequent cloud expansion are ex- 
pected to have no influence on the values of Nmoi and 
No/Njnoi- This is because, by the further ramp, the 
molecules shrink following the field adiabatically. They 
then become tightly bound and weakly interacting among 
themselves, while their local counting is unaffected. Un- 
der these conditions, the condensate fraction, too, should 
not be modified by the subsequent expansion as it is the 
case for an ordinary dilute Bose gas, even though the 
expansion may affect the details of the density profiles. 

In Fig. 4 the condensate fraction No/Nmoi, obtained 
from our theoretical expressions at the lowest tempera- 
ture T/Tc = 0.05, is plotted vs (kpaF)'^ on the BCS 
side of the crossover. The data from Refs. and are 
also reported in the figure. The agreement between the 
overall trends of the theoretical and experimental curves 
appears satisfactory, although quantitative discrepancies 
result between the two sets of curves. They might be 
due to an overestimate of the TF component of the fits 
in Ref. for the reasons discussed in Ref. j9|, and to 
a possible underestimate of the condensate component 
m Ref. I"! due to a preferential loss of molecules in the 
condensate itself. We have, finally, verified that a linear 
dependence occurs between Nmoi/ {Np l'^) and No/Nmoi 
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FIG. 4: Condensate fraction No/Nmoi vs {kpap)^^ on the 
BCS side of the crossover at T /Tc — 0.05 and for the two 
"final" couphngs 0.40 (full line) and 1.50 (dashed line). The 
data from the experiments of Refs. Q (dots) and (squares) 
are also reported. The inset shows Nmoi / (Np /2) vs A'o /Nmoi 
for the same temperature and "final" couplings. 



(inset of Fig. 4) . A similar linear dependence is also re- 
ported in Fig. 4(b) of Ref. although with a different 
definition of Nmoi that includes also molecular states not 
directly detected in the experiment. 
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